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Abstract. A Poincare covariant Faddeev equation, which describes baryons 
as composites of confined-quarks and -nonpointlike-diquarks, is solved to ob- 
tain masses and Faddeev amplitudes for the nucleon and A. The amplitudes 
• are a component of a nucleon-photon vertex that automatically fulfills the 

! Ward-Takahashi identity for on-shell nucleons. These elements are sufficient 

I for the calculation of a quark core contribution to the nucleons' electromag- 

■ netic form factors. An accurate description of the static properties is not 

possible with the core alone but the error is uniformly reduced by the incor- 
poration of meson-loop contributions. Such contributions to form factors are 
noticeable for ^ 2 GeV^ but vanish with increasing momentum transfer. 
Hence, larger experiments probe the quark core. The calculated behaviour 
^ , of G^(Q2)/G^(g2) on Q2 g [2,6] GeV^ agrees with that inferred from po- 

^ I larisation transfer data. Moreover, \/ Q^i*2(Q^)/-^i(Q^) ~ constant on this 

Ch ■ domain. These outcomes result from correlations in the proton's amplitude. 
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1 Introduction 



Modern, high-luminosity experimental facilities that employ large momentum 
transfer reactions are providing remarkable and intriguing new information on 
nucleon structure |l|2j . For an example one need only look so far as the dis- 
crepancy between the ratio of electromagnetic proton form factors extracted via 
Rosenbluth separation and that inferred from polarisation transfer |3I4I5I6I7| . 
This discrepancy is marked for > 2 GeV^ and grows with increasing . At 
such values of momentum transfer, > Af^, where M is the nucleon's mass, a 
veracious understanding of these and other contemporary data require a Poincare 
covariant description of the nucleon. 

This is apparent in applications of relativistic quantum mechanics; e.g., 
Refs. |8l9llUlllll2j . A diff'erent tack follows the formulation of a Poincare co- 
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variant Faddeev equation jl,S|14j . Its foundation is understood through the ob- 
servation that the same interaction which describes colour-singlet mesons also 
generates quark-quark (diquark) correlations in the colour-3 (antitriplet) chan- 
nel ^ni- While diquarks do not survive as asymptotic states pTC]; i.e., they do 
not appear in the strong interaction spectrum, the attraction between quarks in 
this channel grounds a picture of baryons in which two quarks are always corre- 
lated as a colour-3 diquark pseudoparticle, and binding is effected by the iterated 
exchange of roles between the bystander and diquark-participant quarks. 

Reference |T2] reported a rudimentary study of this Faddeev equation and 
subsequently more sophisticated analyses have appeared; e.g., Refs. P^S, 19 20 2T] . 
It has become apparent that the dominant correlations for ground state octet 
and decuplet baryons are scalar and axial-vector diquarks, primarily because the 
associated mass-scales are smaller than the masses of these baryons j22l23j and 
the positive parity of the correlations matches that of the baryons. Both scalar 
and axial- vector diquarks provide attraction in the Faddeev equation; e.g., a 
scalar diquark alone provides for a bound octet baryon and including axial- vector 
correlations reduces that baryon's mass. 

With the retention of axial- vector diquark correlations a quantitative descrip- 
tion of baryon properties is attainable ^Hl- However, that possibility necessitates 
the incorporation of pseudoscalar meson loop contributions because a credible 
description of baryon properties is otherwise problematic. Such effects contribute 
materially: to baryon masses (20^24^25) : and charge and magnetic radii, and mag- 
netic moments |2Hp27j . 

In Sec. 121 we recapitulate on those aspects of the Poincare covariant Faddeev 
equations for the nucleon and A that are important to our analysis, and report 
the solutions. Section describes the formulation of a Ward-Takahashi identity 
preserving current that is appropriate to a nucleon represented by a solution of 
the Faddeev equation. With the necessary elements thus specified, the nucleons' 
electromagnetic form factors are presented and discussed in Sec. H Section El IS 
an epilogue. 

2 Covariant Faddeev Equation 

The properties of light pseudoscalar and vector mesons are described well by 
a rainbow-ladder truncation of QCD's Dyson- Schwinger equations (DSEs) j28j . 
and the calculation of baryon properties using the solution of a Poincare covariant 
Faddeev equation is a desirable extension of that approach. For quarks in the 
fundamental representation of colour-5'[/(3): 

3c 3c 3c = (3c e 6c) 3c = Ic e e 8c © 10c , (1) 

and hence any two quarks in a colour-singlet three-quark bound state must con- 
stitute a relative colour-antitriplet. This enables the derivation of a Faddeev 
equation for the bound state contribution to the three quark scattering ker- 
nel because the same kernel that describes mesons so well is also attractive for 
quark-quark scattering in the colour-3 channel. 

In this truncation of the three-body problem the interactions between two 
selected quarks are added to yield a quark-quark scattering matrix, which is then 
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approximated as a sum over all possible diquark pseudoparticle terms irZJ : Dirac- 
scalar + -axial- vector +[...]. The Faddeev equation thus obtained describes the 
baryon as a composite of a dressed-quark and nonpointlike diquark with an 
iterated exchange of roles between the bystander and diquark-participant quarks. 
The baryon is consequently represented by a Faddeev amplitude: 

^ = ^1+^2 + "^?,^ (2) 

where the subscript identifies the bystander quark and, e.g., ^1^2 are obtained 
from by a correlated, cyclic permutation of all the quark labels. The Faddeev 
amplitude thus obtained does not contain explicit three-body correlations. 

2.1 Ansdtze for the Nucleon and A 

We employ the simplest realistic representation of the Faddeev amplitudes for 
the nucleon and A. The spin- and isospin-1/2 nucleon is a sum of scalar and 
axial-vector diquark correlations: 

MPuC^un)=Mt +M^'', (3) 

with (pi, ai,Ti) the momentum, spin and isospin labels of the quarks constituting 
the bound state, and P = pi + P2 + P3 the system's total momentum. NB. We 
assume isospin symmetry of the strong interaction throughout; i.e., the u- and d- 
quarks are indistinguishable but for their electric charge. Since it is not possible 
to combine an isospin-0 diquark with an isospin-1/2 quark to obtain isospin-3/2, 
the spin- and isospin-3/2 A contains only an axial- vector diquark component 

n^^{pi,ai,n)=Vl\ (4) 

The scalar diquark piece in Eq. Q is 

^ftiP^,c^^,n) = [r^^{lp[,2];K)]zz^'''iK)[sii-,PMP)]z, (5) 

where: ^ the spinor satisfies 

(^7 • P + M) u{P) = = u{P) {i-f-P + M), (6) 

with M the mass obtained by solving the Faddeev equation, and it is also a 
spinor in isospin space with = col(l,0) for the proton and ip- = col(0, 1) 
for the neutron; K = pi + p2 =: P{i2}, P[i2] = Pi - P2, i ■= {-P{12} + 2p3)/3; 
A'^^ is a pseudoparticle propagator for the scalar diquark formed from quarks 
1 and 2, and F^^ is a Bethe-Salpeter-like amplitude describing their relative 
momentum correlation; and S, a 4 x 4 Dirac matrix, describes the relative quark- 
diquark momentum correlation. (5, -T^^ and A^^ are discussed in Sec. 12. 21 ) The 
colour antisymmetry of iZ'a is implicit in P"^ , with the Levi-Civita tensor, eciC2C3) 
expressed via the antisymmetric Gell-Mann matrices; viz., defining 



^The metric we employ is described in [Appendix A:| Euclidean Conventions. 
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then ecic2C3 = {H<'3)cic2- [See Eqs.(jnSl), 

The axial-vector component in Eq. Q is 

where the symmetric isospin-triplet matrices are 

t+ = ^(r° + r') , t° = r\ t- = i^(r° - r') , (9) 

with {T^)ij = dij and t^'^ the usual Pauli matrices, and the other elements in 
Eq. ((SI) are straightforward generalisations of those in Eq. ((SJ • 

The general form of the Faddeev amplitude for the spin- and isospin-3/2 A is 
complicated. However, isospin symmetry means one can focus on the A~^^ with 
it's simple flavour structure, because all the charge states are degenerate, and 
consider 

where Up{P) is a Rarita-Schwinger spinor, Eq. (|A.10|) . 

The general forms of the matrices S{l;P), Al{l;P) and T>up{i;P), which 
describe the momentum space correlation between the quark and diquark in the 
nucleon and the A, respectively, are described in Ref. llQj . The requirement that 
S{l;P) represent a positive energy nucleon; namely, that it be an eigenfunction 
of A+{P), Eq. ()A.7|) . entails 

S{£;P)=si{£;P)lD+(ii-i-i-PlD) S2{t,P), (11) 



where (/D)rs = ^rsi ^ = 1, = —1. In the nucleon rest frame, si^2 describe, 
respectively, the upper, lower component of the bound-state nucleon's spinor. 
Placing the same constraint on the axial- vector component, one has 

6 

^t(^; P) = Y. ^nik P) 75 A^^{1- P) , i = +, 0, - , (12) 

n=l 



where (f^ = I + I ■ P P^, -/^ = -/^ + P P, 



Al = ^-i^P,, Al = -iP,, Al = ^-i^i^, 



(13) 



Finally, requiring also that T>^p{i; P) be an eigenfunction of yl+(P), one obtains 
V,p{e; P) = cS^(Z; P) 5,p + j.A^il; P) , (14) 

with S'^ and given by obvious analogues of Eqs. (fTT|) and (fT^ . respectively. 
One can now write the Faddeev equation satisfied by 'F^ as 



S{k;P)u{P) 
A'^{k;P)u{P) 



M{k,e;P) 



(2vr) 



S{e;P)u{P) 
Aiie;P)uiP) J 



(15) 
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where one factor of "2" appears because ^3 is coupled symmetrically to and ^2, 
and the necessary colour contraction has been evaluated: {H"-)}yc{H"-)cb' = —2 8^1,1. 
The kernel in Eq. is 



with 



0+1 



M{k,e; P) 



MoQ {Mm)i 



(16) 



kq)r' 



\0+, 



where:2 = £ + P/3, kg = k + P/3, 



-qq 



i^q kgq/2; kgg) S {I g) A \^ gg J , 

(17) 

' + 2P/3, kgg = -k + 2P/3 and the 



superscript "T" denotes matrix transpose; and 

{M0l)i=t'rl\kg-£gg/2-£gg) 

XS'^{igg — kg) P^ [£g — k gg / 2] —kgg) S{ig) A 

iM,0t=r^\kg- egg/2; (gg) 
=,1 + 



XS {igg kg)t^r (ig kgq/2; kgg)S{ig)A {^qq) 



{Ml^)t=t^rl\kg-tgg/2;£gg) 



p v^q --qq 

- kg) t* P^ {£g 



kqq/'^'l kgg) S (ig) A]^i^ (igg) 



(18) 
(19) 
(20) 



It is illuminating to note that u(P) in Eq. 1)151) is a normalised average of <p± 
so that, e.g., the proton equation is obtained by projection on the left with 93!^. 
To illustrate this we note that Eq. ()18() generates an isospin coupling between 
u{P)ip+ on the left-hand-side (l.h.s.) of Eq. (|T^ and, on the r.h.s.. 



V2Atu{P)^- -A%iP), 



(21) 



/ip ■^u "-K^ IV+ ■ 

This is just the Clebsch-Gordon coupling of isospin-1© isospin-^ to total isospin- 
^ and means that the scalar diquark amplitude in the proton, {ud)Q+ u, is coupled 
to itself and the linear combination: 

V2 {uu)i+ d — {ud)i+ u . 

Similar statements are obviously true of the spin couplings. 
The Zi's Faddeev equation is 



(22) 



Vxp{k;P) Up{P) 



I 



_dH_ 

(2^ 



Mt{k,l;P)V^^{t,P)u^{P) 



(23) 



with 



M 



* (kq ^qq/'^'Aqq) S {igg kg) t P^ {£g kgg/2; kgg) S (ig) A^J^igg) . 

(24) 



This choice is explained by Eq. 



and the discussion thereabout. 
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2.2 Completing the Faddeev Equation Kernels 

To complete the Faddeev equations, Eqs. (|T5|) & (|23|) . one must specify the 
dressed-quark propagator, the diquark Bethe-Salpeter amphtudes and the di- 
quark propagators that appear in the kernels. 

2.2.1 Dressed-quark propagator 

The dressed-quark propagator can be obtained from QCD's gap equation and 
the general form of the solution is 

S{p) = -i-fpavip^) + as{p^) = l/[n-pA{p'')+B{p^)]. (25) 

It is a longstanding prediction of DSE studies in QCD that the wave function 
renormalisation and dressed-quark mass: 

Z(p2) = 1/A(p2) , M(p2) = i?(p2)/A(p2) , (26) 

respectively, receive strong momentum-dependent corrections at infrared mo- 
menta (29'.:3UI31j : Z(jp) is suppressed and M{p'^) enhanced. The prediction was 
confirmed in recent simulations |S21 of lattice-regularised quenched-QCD,'^ and 
the conditions under which pointwise agreement between DSE results and lattice 
simulations may be obtained have been explored [34f35T36j . The enhancement 
of MijP') is central to the appearance of a constituent-quark mass-scale and an 
existential prerequisite for Goldstone modes. The mass function evolves with in- 
creasing to reproduce the asymptotic behaviour familiar from perturbative 
analyses, and that behaviour is unambiguously evident for > 10 GeV^ |37|38j . 

The importance of this infrared dressing has long been emphasised in DSE 
studies of hadron phenomena and, while numerical solutions of the quark 
DSE are now readily obtained, the utility of an algebraic form for S{p) when 
calculations require the evaluation of numerous multidimensional integrals is 
self-evident. An efficacious parametrisation of S'(p), which exhibits the features 
described above, has been used extensively in hadron studies 128140141] . It is 
expressed via 

as{x) = 27nT{2{x + fi?)) + T{hix)T{h^x)[hQ + h2T{ex)\ , (27) 
av{x) = [1 - ^(2(x + m^))] , (28) 

with X = /}?, fh = m/X, 

Hx) = , (29) 

X 

as{x) = Xasip^) and a'y(x) = }?ay{p'^). The mass-scale, A = 0.566 GeV, and 
parameter values^ 

rn bp hi 62 63 

0.00897 0.131 2.90 0.603 0.185 ' ^ ' 

^For three light flavours the effects of unquenching appear to be small |33|. 

*e = 10"'' in Eq. 1271 acts only to decouple the large- and intermediate-p^ domains. 
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were fixed in a least-squares fit to light-meson observables |42j . The dimensionless 
u = d current-quark mass in Eq. (|5n)) corresponds to 

m„,rf = 5.1MeV. (31) 

The parametrisation yields a Euclidean constituent-quark mass 

M„^, = 0.33GeV, (32) 

defined as the solution of = M^(p^) [3^, whose magnitude is typical of that 
employed in constituent-quark models [H]. In Ref. [l2j, rus = 25 mu^d and = 
0.49 GeV. It is generally true that — M^^ > rhg — rhu^d} where fn denotes 
the renomalisation point independent current-quark mass. The constituent-quark 
mass is an expression of dynamical chiral symmetry breaking, as is the vacuum 
quark condensate^ (^QCD = 0.2 GeV) 

-ml-'^'->?l.^y-^-m2,0.^?. (33) 

Motivated by model DSE studies j4HI47j . Eqs. (l27l) (|28j) express the dressed- 
quark propagator as an entire function. Hence S{p) does not have a Lehmann 
representation, which is a sufficient condition for confinement.^ Employing an en- 
tire function, whose form is only constrained through the calculation of spacelike 
observables, can lead to model artefacts when it is employed directly to calculate 
observables involving large timelike momenta [SII ■ An improved parametrisation 
is therefore being sought. Nevertheless, difficulties are not encountered for mod- 
erate timelike momenta, and on the domain of the complex plane explored in 
the present calculation the integral support provided by an equally effective al- 
ternative cannot differ significantly from that of this parametrisation. 

2.2.2 Diquark Bethe-Salpeter amplitudes 

The rainbow-ladder DSE truncation yields asymptotic diquark states in the 
strong interaction spectrum. Such states are not observed and their appearance 
is an artefact of the truncation. Higher order terms in the quark-quark scattering 
kernel, whose analogue in the quark-antiquark channel do not much affect the 
properties of vector and flavour non-singlet pseudoscalar mesons, ensure that 
QCD's quark-quark scattering matrix does not exhibit singularities which corre- 
spond to asymptotic diquark states ^S]- Nevertheless, studies with kernels that 
do not produce diquark bound states, do support a physical interpretation of 
the masses, 'm(qq)jp , obtained using the rainbow-ladder truncation: the quan- 
tity \qq) jp = ^/''^{qq)jp ^lay be interpreted as a range over which the diquark 
correlation can persist inside a baryon. These observations motivate the Ansatz 

^The condensate is calculated directly from its gauge invariant definition |45| after making 
allowance for the fact that Eqs. (I27II & 12811 yield a chiral-limit quark mass function with 
anomalous dimension 7™ — 1. This omission of the additional ln(p^/ylQQD)-suppression that 
is characteristic of QCD is merely a practical simplification. 

®It is a sufficient condition for confinement because of the associated violation of reflection 
positivity. This notion may be traced from Refs. |46I47I48I49| and is reviewed in Refs. 131140141"] . 
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for the quark-quark scattering matrix that is employed in deriving the Faddeev 
equation: 

[M,,ik, q; i^)]*^ = Yl ^'"(^) ^^"('?; ^) • (34) 

JP=0+,1+,... 

One practical means of specifying the F"^^ in Eq. H34|) is to employ the so- 
lutions of a rainbow-ladder quark-quark Bethe-Salpeter equation (BSE). Using 
the properties of the Gell-Mann matrices one finds easily that := F'^ Ct 
satisfies exactly the same equation as the colour-singlet meson but for a 
halving of the coupling strength p^. This makes clear that the interaction in 
the 3c (qq) channel is strong and attractive.'' Moreover, it follows as a feature of 
the rainbow-ladder truncation that, independent of the specific form of a model's 
interaction, the calculated masses satisfy 

This is a useful guide for all but scalar diquark correlations because the partnered 
mesons in that case are pseudoscalars, whose ground state masses are constrained 
to be small by Goldstone's theorem and which therefore provide a weak lower 
bound. For the correlations relevant herein, models typically give masses (in 
GeV) fim^ : 

^(ud)g+ = 0.74 - 0.82 , m(„„)^^ = m(„rf)^_^ = rrndd)^+ = 0.95 - 1.02 . (36) 

Such values are confirmed by results obtained in simulations of quenched lattice- 
QCD [HJ. Reference |52j also evaluates charge radii for the scalar diquarks: 
r{ud)g+ ~ 1-1 ^'tt; r(^ds)g+ ~ rK+; and r(„^)^^ « 1.3-1.4r^'+. The Bethe-Salpeter 
amphtude is canonically normalised via: 



^ n{K,Q) 



I* 



J 



p 



(37) 
Q=K 



n{K,Q) = trl^r{q;-K)S{q + Q/2)Fiq;K)S^{-q + Q/2). (38) 

A solution of the BSE equation requires a simultaneous solution of the quark- 
DSE 23 . However, since we have already chosen to simplify the calculations by 
parametrising S{p), we follow Ref. j^HI and also employ that expedient with F"' , 
using the following one-parameter forms: 

F^'-{k;K) = -l^H-Ci^^ir^HkV^l+), (39) 

t'Fl\k-K) = ^H-^^^Ct^HkV^l+), (40) 

with the normalisation, J\f'^^, fixed by Eq. ()37() . These Ansdtze retain only that 
single Dirac-amplitude which would represent a point particle with the given 
quantum numbers in a local Lagrangian density: they are usually the dominant 
amplitudes in a solution of the rainbow-ladder BSE for the lowest mass 
diquarks and mesons USESEH- 

^The same analysis shows the interaction to be strong and repulsive in the 6c (qq) channel. 
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2.2.3 Diquark propagators 

Solving for the quark-quark scattering matrix using the rainbow-ladder trun- 
cation yields free particle propagators for A-^ in Eq. (|5H) . As already noted, 
however, higher order contributions remedy that defect, eliminating asymptotic 
diquark states from the spectrum. It is apparent in Ref. ^Hl that the attendant 
modification of A can be modelled efficiently by simple functions that are 
free-particle-like at spacelike momenta but pole-free on the timelike axis. Hence 
we employ 

A'^{K) = ^^(i^V^o^)> (41) 

4:w = (42) 

\ ''''1+ / ''''1+ 

where the two parameters rrijp are diquark pseudoparticle masses and ujjp are 
widths characterising F"^ . Herein we require additionally that 

which is a normalisation that accentuates the free-particle-like propagation char- 
acteristics of the diquarks within the hadron. 

2.3 Nucleon and A Masses 

All elements of the Faddeev equations, Eqs. (P3|l & ((^ . are now completely 
specified. We solve the equations via the method described in Ref. j^S]. The 
masses of the scalar and axial-vector diquarks are the only variable parameters. 
The axial-vector mass is chosen so as to obtain a desired mass for the A, and 
the scalar mass is subsequently set by requiring a particular nucleon mass. 

Two parameter sets are presented in Tabled We obtained Set A by requiring 
a precise fit to the experimental nucleon and A masses. It has long been known 
that this is possible; e.g., Ref. ^\ reports octet and decuplet baryon masses in 
which the rms deviation between the calculated mass and experiment is only 
2%. However, it is also known that such an outcome is undesirable because, 

Table 1. Mass-scale parameters (in GeV) for the scalar and axial- vector diquark correlations, 
fixed by fitting nucleon and A masses: for Set A, a fit to the actual masses was required; 
whereas for Set B the fitted mass was offset to allow for "pion cloud" contributions |2U|. We 
also list LUjp = —r-rrijp, which is the width-parameter in the {qq)jp Bethe-Salpeter amplitude, 
Eqs. 1391 & (1401 : its inverse is a gauge of the diquark's matter radius. 



set 


Mn 


Ma 


mo+ 




UJo+ 




A 


0.94 


1.23 


0.63 


0.84 


0.44=1/(0.45 fm) 


0.59=1/(0.33 fm) 


B 


1.18 


1.33 


0.79 


0.89 


0.56=1/(0.35 fm) 


0.63=1/(0.31 fm) 



1 2 

9 JP 



(43) 
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e.g., studies using the cloudy bag model [M! indicate that the nucleon's mass 
is reduced by as much as SMjsj = —300 to — 400MeV through pion self-energy 
corrections . Furthermore, a perturbative study, using the Faddeev equation, 
of the mass shift induced by pion exchange between the quark and diquark 
constituents of the nucleon obtains (JM^v = —150 to — SOOMeV jUj. We are thus 
led to Set B, which was obtained by fitting to nucleon and A masses that are 
inflated so as to allow for the additional attractive contribution from the pion 
cloud HO]. 

It is apparent in Table ^ that a baryon's mass increases with increasing 
diquark mass, and the fitted diquark mass-scales are commensurate with the 
anticipated values, cf. Eq. (|36]) . with Set B in better accord. If coupling to the 
axial- vector diquark channel is omitted from Eq. (|15|1. then Mff^^ = 1.15 GeV 
and M^^^ = 1.46 GeV. It is thus clear that axial- vector diquark correlations 
provide significant attraction in the nucleon. Of course, using our Faddeev equa- 
tion, the A does not exist without axial-vector correlations. In Set B the amount 
of attraction provided by axial-vector correlations must be matched by that pro- 
vided by the pion cloud. This highlights the constructive interference between the 
contribution of these two effects to a baryons' mass. It is related and noteworthy 
that — mo+ is only a reasonable approximation to — = 0.29 GeV 
when pion cloud effects are ignored: Set A, mi+ — mQ+ = 0.21 GeV cf. Set B, 
mi+ — mQ+ = 0.10 GeV. Plainly, understanding the N-A mass splitting requires 
more than merely reckoning the mass-scales of constituent degrees of freedom. 

3 Nucleon-Photon Vertex 

The nucleon's electromagnetic current is 

J^(P',P) = ieu{P')A^{q,P)u{P), (44) 
= ieu{P') {l^,Fi{Q^) + -i-^^^Q^F2(Q2)^ u{P) , (45) 

where P {P') is the momentum of the incoming (outgoing) nucleon, Q = P' — P, 
and -Fi and F2 are, respectively, the Dirac and Pauli form factors. They are the 
primary calculated quantities, from which one obtains the nucleon's electric and 
magnetic form factors 

GEiQ') = Fi(Q') - ^UQ^) . Gm{Q'') = F,{Q^) + F2(Q2) . (46) 

In Eq. (|44|) . is the nucleon-photon vertex, which we construct following 
the systematic procedure of Ref. This approach has the merit of automat- 
ically providing a conserved current for on-shell nucleons described by Faddeev 
amplitudes of the type we have calculated. Moreover, the canonical normalisa- 
tion condition for the nucleons' Faddeev amplitude is equivalent to requiring 
FiiQ'^ = 0) = 1 for the proton. The vertex has six terms, which are depicted in 



Fig. H and expressed explicitly in Appendix B: Nucleon-Photon Vertex. Here we 



describe the key elements in the construction. 
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Figure 1. Nucleon-photon vertex which ensures a conserved current for on-shell nucleons 
described by the Faddeev ampUtudes, i'i,f, calculated in Sec.|5] The single line represents S{p), 
the dressed-quark propagator, Sec. l2.2.Il and the double line, the diquark propagator, Sec. 12.2.^ 
r is the diquark Bethe-Salpeter amplitude, Sec. 12.2.21 and the remaining vertices are described 
in Secs. l3.lh3.5l the top-left image is diagram 1; the top-right, diagram 2; and so on, with the 
bottom-right image, diagram 6. 



3.1 Diagram 1 

This represents the photon couphng directly to the bystander quark and is ob- 
tained exphcitly from Eqs. ()B.2|) & HB.4|) . It is a necessary condition for current 
conservation that the quark-photon vertex satisfy the Ward-Takahashi identity: 

Qt^ir^iiuh) = S-\ii) - S-\£2) , (47) 

where Q = ii — £2 is the photon momentum flowing into the vertex. Since 
the quark is dressed, Sec. l2.2.ll the vertex is not bare; i.e., /)j(£i,£2) 7^ 7^1- It 
can be obtained by solving an inhomogeneous Bethe-Salpeter equation, which 
was the procedure adopted in the DSE calculation that successfully predicted 
the electromagnetic pion form factor However, since we have parametrised 
S{p), we follow Ref. and write [00] 

i^^{h,i2) = iSA{ilf2)l^. + 2k^ [ij-k^AAiiliD + ABiilil)]; (48) 

with k = {£i+ i2)/2, Q = {£i- I2) and 

SF{£lJl) = \[F{il) + F{il)], AHil4) = ^^^^^, (49) 
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where F = A,B; viz., the scalar functions in Eq. (|25|) . It is critical that in 
Eq. H48() satisfies Eq. (|17|) and very useful that it is completely determined by 
the dressed-quark propagator. Following Ref. [30]) this Ansatz has been used 
fruitfully in many hadronic applications. Its primary defect is the omission of 
pion cloud contributions; e.g., Ref. jSB) but since one of our goals is to draw 
attention to consequences of that omission, this fault is herein a virtue. 

3.2 Diagram 2 

This figure depicts the photon coupling directly to a diquark correlation, for 
which the explicit expression is obtained from Eqs. (|B.2|) & (|B.5|) . In the case of 
a scalar correlation, the general form of the diquark-photon vertex is 

(4, ^2) = 2 f+{k\ k ■ Q, Q2) + /_(fe2, k-Q,Q^), (50) 
and it must satisfy the Ward-Takahashi identity: 

Q^rf (£1,^2) = n^^iij) - n'^iel) , n'^if) = {A-'^if)}-'. (51) 

The evaluation of scalar diquark elastic electromagnetic form factors in Ref. 
is a first step toward calculating this vertex. However, in providing only an 
on-shell component, it is insufficient for our requirements. We therefore adapt 
Eq. (gHl) to this case and write 

F;iUii,h)=k^A^,+ (ele',), (52) 

which is the minimal Ansatz that: satisfies Eq. H51|): is completely determined by 
quantities introduced already; and is free of kinematic singularities. It implements 
/_ = 0, which is a requirement for elastic form factors, and guarantees a valid 
normalisation of electric charge; viz., 

ihn {£', i) = 2e,^ (f) ' =° 2 , (53) 

owing to Eq. (|43|) . NB. We have factored the fractional diquark charge, which 
therefore appears subsequently in our calculations as a simple multiplicative 
factor. 

For the case in which the struck diquark correlation is axial-vector and the 
scattering is elastic, the vertex assumes the form _62_ : ^ 

3 

</?(^i,^2) = -^r^^(^i,^2), (54) 

i=l 

With (r,^(£) = <^a/3-^a^/3/^') 

^i'i(^i'^2) = {h+i2)^.F^xih)Fxf,{e2)F^iil4), (55) 
rj^U^uh) = [T^a{h)T^p{h)ei, + T^f,{e2)T^p{h)hp]F2{£l£h, (56) 

^ia/3(^l'^2) = -:7^{il+£2),.T^p{h)i2pTpxii2)ilxFs{ilel). (57) 

*If the scattering is inelastic the general form of the vertex involves eight scalar functions |63|. 
Absent further constraints and input, we ignore the additional structure in this Ansatz. 
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This vertex satisfies: 

hc.rllp{ha2) = = rllp{i^a2)i2p , (58) 

which is a general requirement of the elastic electromagnetic vertex of axial- 
vector bound states and guarantees that the interaction does not induce a pseu- 
doscalar component in the axial-vector correlation. We note that the electric, 
magnetic and quadrupole form factors of an axial-vector bound state are ex- 
pressed 

G'/iQ') = Fi + f ri+ G'^iQ') , n+ = (59) 

Gjv,(Q') = -F2(Q') , (60) 
G'q{Q') = FiiQ^) + F^iQ') + (1 + ri+) FsiQ') . (61) 

Extant knowledge of the form factors in Eqs. (|54|) " (|57j) is limited and thus 
one has little information about even this rudimentary vertex model. Hence, we 
employ the following Ansdtze: 

F,{£lil) = A^,^{£l£l), (62) 
F2iil4) = -Fi + (l-ri+)(ri+Fi + l-/ii+)d(ri+) (63) 
mlJl) = -(xi+(l-ri+)(i(ri+)+Fi+F2)d(ri+), (64) 

with d{x) = 1/(1 + x)^. This construction ensures a valid electric charge normal- 
isation for the axial- vector correlation; viz., 

lun r^i:^(/, i) = T^pii) A 7ji+ (^2) i'-o ^^^^^^ 2 , (65) 

owing to Eq. ()43() . and current conservation 

nl+ 

The diquark's static electromagnetic properties follow: 

(0) = 1 , (0) = , G'^ (0) = . (67) 

For a pointlike axial-vector: = 2; and Xi+ = 1) which corresponds to an 
oblate charge distribution. In addition, Eqs. ((5Hl - H57|) with Eqs. H62() -(|M |) realise 
the constraints of Ref. [S^; namely, independent of the values of fii+ & Xi+i t^is 
form factors assume the ratios 

Gl\Q') : G^(Q2) : G^^ {Q^) ^'=°^ (1 - |n+) : 2 : -1 . (68) 
3.3 Diagram 3 

This image depicts a photon coupling to the quark that is exchanged as one 
diquark breaks up and another is formed. While this is the first two-loop diagram 



hm Q,rl^p{h,£2) = 0. (66) 
£2— 
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we have described, no new elements appear in its specification: the dressed-quark- 
photon vertex was discussed in Sec. 13. II The exphcit expression for this diagram's 
contribution to the nucleons' form factors is obtained from Eqs. I|B.2|) & ()B.7|) . 

It is noteworthy that the process of quark exchange provides the attraction 
necessary in the Faddeev equation to bind the nucleon. It also guarantees that 
the Faddeev amplitude has the correct antisymmetry under the exchange of 
any two dressed-quarks. This key feature is absent in models with elementary 
(noncomposite) diquarks. 

3.4 Diagram 4 

This differs from Diagram 2 in expressing the contribution to the nucleons' form 
factors owing to an electromagnetically induced transition between scalar and 
axial-vector diquarks: the explicit expression is described in connection with 
Eq. I|B.5|) . The transition vertex is a rank-2 pseudotensor, kindred to the matrix 
element describing the p^*Tr^ transition and can therefore be expressed 

rZi^uh) = -rTsi^uh) = j^T{iui2)e^^pxhpi2X, (69) 

where 7, a are, respectively, the vector indices of the photon and axial-vector 
diquark. For simplicity we proceed under the assumption that 

Tih,i2) = Kr; (70) 

viz., a constant, for which a typical value is |66j : 

~ 2 . (71) 

In the nucleons' rest frame, a conspicuous piece of the Faddeev amplitude 
that describes an axial- vector diquark inside the bound state can be characterised 
as containing a bystander quark whose spin is antiparallel to that of the nucleon, 
with the axial-vector diquark's parallel. The interaction pictured in this diagram 
does not affect the bystander quark but the transformation of an axial-vector 
diquark into a scalar effects a flip of the quark spin within the correlation. After 
this transformation, the spin of the nucleon must be formed by summing the 
spin of the bystander quark, which is still aligned antiparallel to that of the 
nucleon, and the orbital angular momentum between that quark and the scalar 
diquark.^ This argument, while not sophisticated, does motivate an expectation 
that Diagram 4 will strongly impact on the nucleons' magnetic form factors. 

3.5 Diagrams 5 &: 6 

These two-loop diagrams are the so-called "seagull" terms, which appear as part- 
ners to Diagram 3 and arise because binding in the nucleons' Faddeev equations 

®A less prominent component of the amplitude has the bystander quark's spin parallel to that 
of the nucleon while the axial- vector diquark's is antiparallel: this © (55) j+ system has 
one unit of angular momentum. That momentum is absent in the ® (99)0+ system. Other 
combinations also contribute via Diagram 3 but all mediated processes inevitably require a 
modification of spin and/or angular momentum. 
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is effected by the exchange of nonpointlike diquark correlations [2H! • The exphcit 
expression for their contribution to the nucleons' form factors is obtained from 
Eqs. (|B.2|) and ()B.8|) . The new elements in these diagrams are the couplings of a 
photon to two dressed-quarks as they either separate from (Diagram 5) or com- 
bine to form (Diagram 6) a diquark correlation. As such they are components 
of the five point Schwinger function which describes the coupling of a photon to 
the quark-quark scattering kernel. This Schwinger function could be calculated, 
as is evident from the recent computation of analogous Schwinger functions rele- 
vant to meson observables |2.7J. However, such a calculation provides valid input 
only when a uniform truncation of the DSEs has been employed to calculate 
each of the elements described hitherto. We must instead employ an algebraic 
parametrisation [HH], which for Diagram 5 reads 



Xf (A;,Q) = eby 



+ ee 



4fc^ - 
Ak-Q-Q^ L 

4A; • Q + Q2 



rJ^'ik - Q/2) - r-^^'ik) 
r^^k + Q/2) - r^\k) 



(72) 



with k the relative momentum between the quarks in the initial diquark, Cby the 
electric charge of the quark which becomes the bystander and eox, the charge of 
the quark that is reabsorbed into the final diquark. Diagram 6 has 



xr. 



- e» 



Ak^ - Q^ 
Ak-Q-Q^' 
Akf, + Q^ 
Ak-Q + Q^ 



r 



{k + Q/2) - r 

J\k - Q/2) - 



r'\k) 



(k) 



(73) 



— tP 

where F 



{£) is the charge-conjugated amplitude, Eq. HA.9|) . Plainly, these terms 
vanish if the diquark correlation is represented by a momentum-independent 
Bethe-Salpeter-like amplitude; i.e., the diquark is pointlike. 

It is naturally possible to use more complicated Ansdtze. However, like 
Eq. H52|) . Eqs. (|72|) & ()73() are simple forms, free of kinematic singularities and 
sufficient to ensure the nucleon-photon vertex satisfies the Ward-Takahashi iden- 
tity when the composite nucleon is obtained from the Faddeev equation. 



4 Nucleon Electromagnetic Form Factors 



4.1 Exegesis 

In order to place the calculation of baryon observables on the same footing as the 
study of mesons, the proficiency evident in Refs. j28!(i7J will need to be applied 
to every line and vertex that appears in Fig.^ While that is feasible, it requires 
a considerable investment of resources. In the meantime, herein we present a 
study whose merits include a capacity to: explore the potential of the Faddeev 
equation truncation of the baryon three-body problem; and elucidate the role of 
additional correlations, such as those associated with pseudoscalar mesons. 

It is worthwhile to epitomise our input before presenting the results. One 
element is the dressed-quark propagator, Sec. 12. 2.11 The form we use both 
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anticipated and expresses the features that are now known to be true It 
carries no free parameters, because its behaviour was fixed in analyses of meson 
observables, and is basic to a description of hght- and heavy-quark mesons that 
is accurate to better than 10% P7| . 

We proposed that the nucleon is at heart composed of a dressed-quark and 
nonpointlike diquark with binding effected by an iterated exchange of roles be- 
tween the bystander and diquark-participant quarks. This picture is realised via 
a Poincare covariant Faddeev equation, Sec. 12. 11 which incorporates scalar and 
axial- vector diquark correlations. There are two parameters, Secs. l2.2?2] fc 12.2^ 
the mass-scales associated with these correlations. They are fixed by fitting to 
specified nucleon and A masses, Sec. l2.3| and thus at this point there are still no 
free parameters with which to influence the nucleons' form factors. 

With the constituents and the bound states' structure defined, only a spec- 
ification of the nucleons' electromagnetic interaction remained. Its formulation 
was guided almost exclusively by a requirement that the nucleon-photon vertex 
satisfy a Ward-Takahashi identity. Since the scalar diquark's electromagnetic 
properties are readily resolved, our result, Fig.^ depends on three parameters 
that are all tied to properties of the axial- vector diquark correlation: & xi+ ) 
respectively, the axial-vector diquarks' magnetic dipole and electric quadrupole 
moments; and Kq-, the strength of electromagnetic axial- vector <-> scalar diquark 
transitions. Hence, with our calculations we exhibit and interpret the dependence 
of the nucleons' form factors on these three parameters, and also on the nucle- 
ons' intrinsic quark structure as expressed in the Poincare covariant Faddeev 
amplitudes. 



4.2 Calculated Results and Discussion 
4.2.1 Static properties and form factors 

The nucleons' charge and magnetic radii are defined respectively as 



r% :=-6^1nGf(s) 



s=0 



(74) 



s=Q 



where N = n,p; and in Table EJ we report charge radii calculated for a range 
of values of the parameters that characterise the axial- vector diquarks' electro- 
magnetic form factors. Sec. 13. 21 centred on the point-particle values of fii+ = 2 
& Xi+ = 1) Eg. (|67jl . and Kq- = 2, Eg. (|7l|) . The radii, particularly that of 
the neutron, are most sensitive to changes in the axial-vector diguarks' elec- 
tric guadrupole moment, Xi+- This is not surprising given that Xi+ is only 
model parameter that speaks directly of the axial- vector diquarks' electric charge 
distribution. The radii's insensitivity to kt, the strength of the scalar <-> axial- 
vector transition, is concordant with the discussion in Sec. 13. 41 With the reference 
values given in Eqs. (|H7|) & (fTT]) . Set A underestimates the proton radius by 30% 
and the magnitude of the neutron radius by 43%, while for Set B these differences 
are 32% and 50%, respectively. 

Table El presents results for the nucleons' magnetic radii. They are insensi- 
tive to the axial-vector diquarks' quadrupole moment but react to the diquarks' 
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Table 2. Charge radii, in fm, calculated using the diquark mass-scale parameters in Table Q 
for a range of axial-vector-diquark-photon vertex parameters, centred on the point-particle 
values of ^]^+ = 2 k, Xi+ = 1, Eq. H67^ . and kt = 2, Eq. H71^ . Columns labelled o give the 
percentage-difference from results obtained with the reference values. r„ := — ^ — {r^) . Values 
inferred from experiment are (58 : rp — 0.847 & r„ = —0.336. 

Set A Set B 





Xi+ 






rp 


rn 


' n 


rp 


rp 


rn 


' n 


1 


1 


2 


0.599 


-1.2 


0.185 


-4.1 


0.596 


0.2 


0.171 


1.2 


2 


1 


2 


0.606 




0.193 




0.595 




0.169 




3 


1 


2 


0.614 


1.3 


0.200 


3.6 


0.593 


-0.3 


0.167 


-1.2 


2 





2 


0.593 


-2.2 


0.179 


-7.3 


0.575 


-3.4 


0.145 


-14.2 


2 


2 


2 


0.620 


2.3 


0.205 


6.2 


0.614 


3.2 


0.191 


13.0 


2 


1 


1 


0.606 


0.0 


0.189 


-2.1 


0.595 


0.0 


0.167 


-1.2 


2 


1 


3 


0.606 


0.0 


0.196 


1.6 


0.595 


0.0 


0.172 


1.8 



Table 3. Magnetic radii, in fm, calculated with the diquark mass-scales in Tableland the 
parameter range described in Table |5| Columns labelled a give the percentage-difference from 
results obtained with the reference values: = 2, Xi+ — li '^t = 2. Values inferred from 
experiment are = 0.836 & < = 0.889. 



Set A Set B 





Xi+ 




rp 


<^ 

rp 


rn 


' n 


rp 




rn 


' n 


1 


1 


2 


0.456 


-2.4 


0.467 


-1.3 


0.442 


-1.6 


0.446 


-0.7 


2 


1 


2 


0.467 




0.473 




0.449 




0.449 




3 


1 


2 


0.477 


2.1 


0.478 


1.1 


0.454 


1.1 


0.454 


1.1 


2 





2 


0.467 


0.0 


0.473 


0.0 


0.449 


0.0 


0.449 


0.0 


2 


2 


2 


0.467 


0.0 


0.473 


0.0 


0.449 


0.0 


0.449 


0.0 


2 


1 


1 


0.470 


0.6 


0.480 


1.3 


0.453 


0.9 


0.459 


2.2 


2 


1 


3 


0.465 


-0.4 


0.472 


-0.2 


0.445 


-0.9 


0.446 


-0.7 



magnetic moment as one would anticipate: increasing in magnitude as in- 
creases. Moreover, consistent with expectation, Sec. 13. 41 these radii also respond 
to changes in Kq-, decreasing as this parameter is increased. With the reference 
values in Eqs. H67() & (fTT]) . both Sets underestimate by approximately 40%. 

TableEllists results for the nucleons' magnetic moments. They, too, are insen- 
sitive to the axial-vector diquarks' quadrupole moment but react to the diquarks' 
magnetic moment, increasing quickly in magnitude as increases. As antici- 
pated in Sec. 13. 41 the nucleons' moments respond strongly to alterations in the 
strength of the scalar <-> axial-vector transition, increasing rapidly as kj- is in- 
creased. Set A, which is fitted to the experimental values of Mn k. M/\, describes 
the nucleons' moments quite well: Kp is 15% too large; and \^n\^ 16% too small. 
On the other hand, Set B, which is fitted to baryon masses that are inflated so as 
to make room for pion cloud effects, overestimates Kp by 47% and by 18%. 
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Table 4. Magnetic moments, in nuclear magnetons, calculated with the diquark mass-scales in 
Tableland the parameter range described in Table |5| Columns labelled a give the percentage- 
difference from results obtained with the reference values: = 2, Xi+ = 1, i^t = 2. Experi- 
mental values are: ~ fj,p — 1 = 1.79 &z. fj,„ — —1.91. 



Set A Set B 



^1+ 


Xi+ 












Kp 




l^nl 




1 


1 


2 


1.79 


-15.3 


1.70 


-5.1 


2.24 


-21.9 


2.00 


-6.2 


2 


1 


2 


2.06 




1.79 




2.63 




2.13 




3 


1 


2 


2.33 


15.4 


1.88 


5.1 


3.02 


21.9 


2.26 


6.1 


2 





2 


2.06 


0.0 


1.79 


0.0 


2.63 


0.0 


2.13 


0.0 


2 


2 


2 


2.06 


0.0 


1.79 


0.0 


2.63 


0.0 


2.13 


0.0 


2 


1 


1 


1.91 


-8.4 


1.64 


-8.4 


2.45 


-10.1 


1.95 


-8.5 


2 


1 


3 


2.21 


8.4 


1.85 


8.3 


2.82 


10.1 


2.31 


8.5 



Nucleon electromagnetic form factors associated with the tabulated values 
of static properties are presented in Figs.[2Hll These figures confirm and aug- 
ment the information in Tables [2H11 Consider, e.g., the electric form factors. 
One observes that the differences between results obtained with Set A and Set B 
generally outweigh those delivered by variations in the parameters characterising 
the axial-vector diquark's electromagnetic properties. The proton's electric form 
factor, in particular, is largely insensitive to these parameters, and it is apparent 
that the nucleon's electric form factor only responds notably to variations in 
Xi+, Fig. 1^1 The nucleons' magnetic form factors exhibit the greatest sensitivity 
to the axial-vector diquarks's electromagnetic properties but in this case, too, 
the differences between Set A and Set B are more significant. For > 4GeV^ 
there is no sensitivity to the diquarks' electromagnetic parameters in any curve. 
This is naturally because our parametrisation expresses the perturbative limit, 
Eq. H68|) . It is thus apparent from these figures that the behaviour of the nu- 
cleons' form factors is primarily determined by the information encoded in the 
Faddeev amplitudes. 

Our results show that the nucleons' electromagnetic properties are sensitive to 
the strength of axial- vector diquark correlations in the bound state and react to 
the electromagnetic properties of these correlations. In all cases the dependence is 
readily understood intuitively. However, taken together our results indicate that 
one cannot readily tune the model's parameters to provide a uniformly good 
account of nucleon properties: something more than dressed-quark and -diquark 
degrees of freedom is required. 

4.2.2 Chiral corrections 

It is appropriate now to examine effects that arise through coupling to pseu- 
doscalar mesons. As with baryon masses, there are two types of contribu- 
tions from meson loops to electromagnetic form factors: regularisation-scheme- 
dependent terms, which are analytic in the neighbourhood of m = 0; and nonana- 
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[GeV^] [Gev"] 

Figure 2. Response of nucleon form factors to variations in the magnetic moment of the axial- 
vector diquark: = 1, 2, 3; with Xi+ = 1, ~ 2. The legend in the top-left panel applies to 
all; the dashed-line marked by "*" is a fit to experimental data (69| and the dashed-line marked 
by in the lower- left panel is the fit to G%{Q^) of Ref. I7U| : and the horizontal lines in the 
right panels mark the experimental value of the nucleon's magnetic moment. 



lytic scheme-independent terms. For the static properties presented in Tables[2]-El 
the leading-order scheme-independent contributions are [7J 



where gA = 1-26, = 0.0924 GeV= 1/(2. 13fm), = fip — Hn- Clearly, the radii 
diverge in the chiral limit, a much touted aspect of chiral corrections. 

While these scheme-independent terms are immutable, at physical values 
of the pseudoscalar meson masses they do not usually provide the dominant 
contribution to observables: that is provided by the regularisation-parameter- 
dependent terms. This is apparent for baryon masses in Ref. [20] and for the 
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[GeV^] [Gev"] 

Figure 3. Response of nucleon form factors to variations in the electric quadrupole moment of 
the axial-vector diquark: Xi+ ~ 0, 1, 2; with = 2, kt = 2. The other features are described 
in the caption of Fig.|5| 

pion charge radius in Ref. 61 . It is particularly important here, as is made plain 
by a consideration of the neutron charge radius. From Eq. 1)75(1 . one obtains 

(-n)]v1°°' = -(0.48fm)2, (78) 

which is more than twice the experimental value. On the other hand, the contri- 
bution from the low energy constants is [7J 

(-n)L'J°°' = +(0-69fm)2, (79) 

which is four-times larger in magnitude than the experimental value and has the 
opposite sign. This emphasises the delicate cancellation that is arranged in chiral 
perturbation theory to fit the neutron's charge radius. In this instance the re- 
maining important piece is the neutron's Foldy term: 3/i„/(2M^) = — (0.35fm)^, 
which is also a fitted quantity. Moreover, for the magnetic radii it is established 
that, at the physical pion mass, the leading chiral limit behaviour is not a good 
approximation [JJ . Additional discussion of issues that arise in formulating a chi- 
ral expansion in the baryon sector and its convergence may be found in Ref. |72j . 
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Figure 4. Response of nucleon form factors to variations in the strength of the electromagnetic 
a:xial- vector ^ scalar diquark transition: kt = 1, 2, 3; with = 2, Xi+ = 1- The other features 
are described in the caption of Fig.|5| 



Since regularisation-parameter-dependent parts of the chiral loops are im- 
portant we follow Ref. j73j and estimate the corrections using modified formulae 
that incorporate a single parameter which mimics the effect of regularising the 
integrals. Thus Eqs. ((75|) & (f7S|) are rewritten 

ir\)Vr' - ^'^^M^), (80) 



il^^JNA = ^ A 2f2 "^TT-arctan — 



(81) 
(82) 



wherein A is a regularisation mass-scale, for which a typical value is ~ 0.4 GeV 
[73]. NB. The loop contributions vanish when the pion mass is much larger than 
the regularisation scale, as required: very massive states must decouple from 
low-energy phenomena. 
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Table 5. Row 1 - static properties calculated with Set B diquark masses, Tabled and — 2, 
Xi+ = 1, Kr = 2; charge radii in fm, with r-„ defined in Table |5| and magnetic moments in 
nuclear magnetons. Row 2 adds the corrections of Eqs. H8()^ - H82|l with A = 0.3 GeV. c in row n, 
is the rms relative-difference between the entries in row n and 3. 





rp 












<? 


q-{qq) core 


0.595 


0.169 


0.449 


0.449 


3.63 


2.13 


0.39 


-|-7r-loop correction 


0.762 


0.506 


0.761 


0.761 


3.05 


1.55 


0.23 


experiment 


0.847 


0.336 


0.836 


0.889 


2.79 


1.91 





We return now to the calculated values of the nucleons' static properties, 
Tables and focus on the Set B results obtained with = 2, Xi+ = 1> 
Kq- = 2. Recall that Set B was chosen to give inflated values of the nucleon and 
A masses in order to make room for chiral corrections, and therefore one may 
consistently apply the corrections in Eas. ((771) . (|Sn|) & H81() to the static prop- 
erties. With A = 0.3 GeV this yields the second row in Table El the regularised 
chiral corrections reduce the rms relative-difference signficantly. This crude anal- 
ysis, complementing Ref. [SO], suggests that a veracious description of baryons 
can be obtained using dressed-quark and -diquark degrees of freedom augmented 
by a sensibly regulated pseudoscalar meson cloud. 



4.2.3 Form factor ratios 

In Fig.Elwe plot the ratio ^ip G^(Q^)/G^^(Q2). The behaviour of the experimen- 
tal data at small is readily understood. In the neighbourhood of = 0, 

and because Vp Vp the ratio varies by less than 10% on < < 0.6 GeV^, if 
the form factors are approximately dipole. This is evidently true of the experi- 
mental data. In our calculation, without chiral corrections. Table |^&: 01 rp > rp. 
Hence the ratio must fall immediately with increasing Q^. Incorporating pion 
loops, we obtain the results in Row 2 of TablejSl which have rp ^ rp. The small 

behaviour of this ratio is thus materially affected by the proton's pion cloud. 

True pseudoscalar mesons are not pointlike and therefore pion cloud contribu- 
tions to form factors diminish in magnitude with increasing Q^. For example, in 
a study of the A transition |7S1, pion cloud contributions to the Ml form 

factor fall from 50% of the total at = q to < 10% for > 2GeV2. Hence, 
the evolution of ^p G^^{Q'^) / G\j{Q'^) on > 2GeV^ is primarily determined 
by the quark core of the proton. This is evident in Fig.jSl which illustrates that, 
on E (1, 5) GeV^, fip G^^{Q'^) / G\j{Q'^) is sensitive to the parameters defining 
the axial-vector-diquark-photon vertex. The response diminishes with increasing 

because our parametrisation expresses the perturbative limit, Eq. (|68j) . 

The behaviour of fipG^j^{Q^)/G^lj{Q^) on > 2 GeV^ is determined either 
by correlations expressed in the Faddeev amplitude, the electromagnetic prop- 
erties of the constituent degrees of freedom, or both. The issue is decided by the 
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Figure 5. Proton form factor ratio: fip G^{Q'^)/G^{Q'^). Calculated results: lower band - Set A 
in Table and upper band - Set B. For both bands, G^{Q^) was calculated using the point- 
particle values: = 2 & Xi+ ~ li Eq. 1671 . and kt = 2, Eq. 17111 : i.e., the reference values 
in Tables HEl Variat ions in the axial- vector diquark parameters used to evaluate G^((5^) have 
little effect on the plotted results. The width of the bands reflects the variation in G\j{Q^) with 
axial- vector diquark parameters and, in both cases, the upper border is obtained with = 3, 
Xi+ ~ 1 and K.T = 2, while the lower has = 1. The data are: squares - Ref. diamonds - 
Ref. 0; and circles - Ref. |71). 




fact that the magnitude and trend of the results are not materially affected by 
the axial- vector-diquarks' electromagnetic parameters. This observation suggests 
strongly that the ratio's evolution is due primarily to spin-isospin correlations in 
the nucleon's Faddeev amplitude. One might question this conclusion, and argue 
instead that the difference between the results depicted for Set A cf. Set B origi- 
nates in the larger quark-core nucleon mass obtained with Set B (Table ^) which 
affects through Eq. We checked and that is not the case. Beginning 

with the Set B results for F\^2 we calculated the electric form factor using the 
Set A nucleon mass and then formed the ratio. The result is very different from 
the internally consistent Set A band in Fig.O e.g., it drops more steeply and lies 
uniformly below, and crosses zero for ~ 3.7 GeV^. 

It is noteworthy that Set B, which anticipates pion cloud effects, is in reason- 
able agreement with both the trend and magnitude of the polarisation transfer 
data |3|4l5j . NB. Neither this nor the Rosenbluth |74j data played any role in 
the preceding analysis or discussion. 
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Q' [GevT 

Figure 6. Proton Pauli/Dirac form factor ratios. The data are as described in Fig.|3 as are 
the bands except that here the upper border is obtained with = 1, Xi+ = 1 a-nd nr ~ 2, 
and the lower with =3. 



We have also examined the proton's Dirac and Pauh form factors in isolation. 
On the domain covered, neither Fi(Q'^) nor i^2(Q^) show any sign they have 
achieved the asymptotic behaviour anticipated from perturbative QCD. 

In Fig.lSJwe depict weighted ratios of these form factors. Our numerical results 
are consistent with 

VQ^^^ ^ constant, 2 < ^^(GeV^) < 6, (84) 

as are the polarisation transfer data. Such behaviour has been argued to indi- 
cate the presence of substantial orbital angular momentum in the proton |8|76j . 
Orbital angular momentum is not a Poincare invariant. However, if absent in 
a particular frame, it will almost inevitably appear in another frame related 
via a Poincare transformation. Nonzero quark orbital angular momentum is the 
necessary outcome of a Poincare covariant description. This is why the covari- 
ant Faddeev amplitude describing a nucleon is a matrix-valued function with 
a rich structure that, in the nucleons' rest frame, corresponds to a relativistic 
wave function with s-wave, p-wave and even d-wave components. (Details can 
be found in Ref. i77j. Sec. 2.4.) The result in Fig.|21is not significantly influenced 
by details of the diquarks' electromagnetic properties. Instead, the behaviour is 
primarily governed by correlations expressed in the proton's Faddeev amplitude 
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Figure 7. Weighted proton Pauli/Dirac form factor ratio, calculated with A = Mjv = 
0.94 GeV. The bands are as described in Fig.|Sl as are the data. 

and, in particular, by the amount of intrinsic quark orbital angular momentum 
j59j ■ NB. This phenomenon is analogous to that observed in connection with the 
pion's electromagnetic form factor. In that instance axial-vector components of 
the pion's Bethe-Salpeter amplitude are responsible for the large behaviour 
of the form factor: they alone ensure Q'^F.,^{Q'^) ~ constant for truly ultraviolet 
momenta (2Sj- These components are required by covariance ,45^ and signal the 
presence of quark orbital angular momentum in the pseudoscalar pion. 

In Fig.EI we plot another weighted ratio of Pauli and Dirac form factors. 
A perturbative QCD analysis |79j that considers effects arising from both the 
proton's leading- and subleading-twist light-cone wave functions, the latter of 
which represents quarks with one unit of orbital angular momentum, suggests 



where 71 is a mass-scale that corresponds to an upper-bound on the domain of 
nonperturbative (soft) momenta. This scaling hypothesis is not predictive unless 
the value of A is known a priori. However, A cannot be computed in perturbation 
theory. 

A scale of this type is not an elemental input to our calculation. It is instead a 
derivative quantity that expresses the net integrated effect of many basic features, 
among which are the mass-scale characterising quark-dressing and that implicit 




constant, ^ A^ , 



(85) 
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in the support of the Faddeev amphtude. Extending our calculation to larger 
is not a problem in principle, and the numerical challenge may readily be met 
by choosing to use something more powerful than a single desktop computer. 
With such information we would be confident of credibly exploring the validity 
of Eq. (|85|) : lacking it, we can only provide a plausible argument. 

One needs an estimate of a reasonable value for A. The nucleon's mass, M^r, 
is one natural mass-scale in our calculation. Other relevant mass-scales are those 
which characterise the electromagnetic size of the nucleon and its constituents. A 
dipole mass-scale for the proton is approximately 0.85 GeV; the dressed-quark- 
photon- vertex is characterised by a monopole mass-scale of 0.8 GeV [SHI; and 
the diquark-photon vertices, by monopole mass-scales ^/Srrijp ~ 1.0- 1.5 GeV, 
Eqs. ((^ & 

As an adjunct one can consider the dressed-quark mass function, defined in 
Eq. ((2ni) and discussed thereabout. A nonzero mass function in the chiral limit 
is an essentially nonperturbative phenomenon. Hence the ratio 



vanishes on the perturbative-Q^ domain. For = 0, on the other hand, cal- 
culations typically yield j38j -Rj^(O) = 0.96; i.e., the mass function's behaviour 
is almost completely nonperturbative. The Q^-evolution of R^{Q'^) can there- 
fore guide in demarcating the nonperturbative domain. Reference |34j provides a 
mass-function that agrees pointwise with quenched-QCD lattice data and gives 
a unique chiral-limit mass function. From these results one finds that R'^Q^ > 
4M^) < 0.5i?2(0)) '^i^-! perturbative effects are dominant in the u-quark mass 
function on £ (4M^, oo). On the other hand, 0.8 < < M^) and hence 

the li-quark mass function is principally nonperturbative on E [0,M^]. 

Together these observations suggest that, in our model and in QCD, a ju- 
dicious estimate of the least-upper-bound on the domain of soft momenta is 
A = Mjsr, and this is the value employed in Fig.d the figure does not provide 
compelling evidence for Eq. H85() on the domain for which information is currently 
available. One can attempt to fit Eq. (|H^ to the calculated results or data, and 
[TH] A = ylfit « 0.3 GeV, provides fair agreement. However, this value is sig- 
nificantly smaller than the natural scales we have identified in the model and 



5 Epilogue 

We recapitulated on a Poincare covariant Faddeev equation that describes 
baryons as composites of confined-quarks and -diquarks, and solved this equation 
to obtain masses and amplitudes for the nucleon and A. Two parameters appear 
in the model Faddeev equation: masses of the scalar and axial-vector diquark 
correlations. They were fixed by fitting stipulated masses of the baryons. We 
interpreted the masses and Faddeev amplitudes thus obtained as representing 
properties of the baryons' "quark core," and argued that this should be aug- 
mented in a consistent fashion by chiral-loop corrections. 



(86) 



i?0(4J = 0.98i?0(0). 
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We explained subsequently the formulation of a nucleon-photon vertex, which 
automatically ensures the vector Ward-Takahashi identity is fulfilled for on-shell 
nucleons described by the calculated Faddeev amplitudes. This guarantees cur- 
rent conservation. The vertex Ansatz involves three parameters. Two of these 
specify electromagnetic properties of axial-vector diquarks and a third measures 
the strength of electromagnetically induced axial-vector- scalar-diquark tran- 
sitions. These quantities are also properties of the nucleons' quark core. 

The elements just described arc sufficient for a calculation of the quark contri- 
bution to the nucleons' electromagnetic form factors. We explored a reasonable 
range of nucleon-photon-vertex parameter values and found that an accurate 
description of the nucleons' static properties was not possible with the core com- 
ponents alone. However, this mismatch with experiment was greatly reduced by 
the inclusion of chiral corrections. 

We calculated ratios of the proton's form factors. On the whole domain of 
nucleon-photon-vertex parameter values explored, the calculated behaviour of 
G^{Q'^)/G^^{Q'^) for > 2GeV^ agrees with that inferred from contemporary 
polarisation transfer data. Moreover, with the same insensitivity to parameters, 
the ratio ■y/Q^F2{Q'^)/Fi{Q'^) constant on G [2, 6] GeV^. Since the parame- 
ters in the nucleon-photon vertex do not influence these outcomes, we judge they 
are manifestations of features intrinsic to the micleon's Faddeev amplitude. In the 
nucleon's rest frame, this amplitude corresponds to a relativistic wave function 
with S-, p- and even d-wave quark orbital angular momentum components. 

In our view baryons can realistically be seen as a dominant Poincare covari- 
ant quark core, augmented by pseudoscalar meson cloud contributions that, e.g., 
make a noticeable contribution to form factors for ^ 2 GeV^ . Meson com- 
positeness ensures that such contributions diminish with increasing Q^. Hence 
experiments at larger serve as an instructive probe of correlations in baryon 
wave functions; i.e., their Faddeev amplitudes. A good understanding of QCD's 
long-range dynamics is required in order to obtain a reliable quark-core wave 
function. This is the box that contemporary experiments are opening. 
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Appendix A: Euclidean Conventions 

In our Euclidean formulation: 

4 

p - q = ^ Piqi; (A.l) 

i=l 

{Tm. li-} = 2 5^,.; 7^ = 7p ; Oy.^ = | [^f^f^] ; tr[767^7i,7p7a] = -4 e^upa , £1234 = 1 - (A.2) 
A positive energy spinor satisfies 

u{P, s){i-fP + M) = 0={i'yP + M) u{P, s) , (A.3) 
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where s = ± is the spin labeL It is normahsed: 

u(P, s) u{P, s) = 2M (A.4) 

and may be expressed exphcitly: 

u{P, s) = VM - i£ { ^a P ^ ) , (A.5) 



with £ = iVP^ + M^ 



M - i£ ' 



(A.6) 



For the free-particle spinor, u(P, s) — u{P, 3)^74. 

The spinor can be used to construct a positive energy projection operator: 



s — ± 

A negative energy spinor satisfies 

v{P, s) (i7 • P - Af) = = (i7 • P - M) v{P, s) , (A.8) 

and possesses properties and satisfies constraints obtained via obvious analogy with u(P, s). 
A charge-conjugated Bethe-Salpeter amplitude is obtained via 

r{k-P)=C''r{-k;PfC, (A.9) 

where "T" denotes a transposing of all matrix indices and C — 7274 is the charge conjugation 
matrix, = — C. 

In describing the A resonance we employ a Rarita-Schwinger spinor to unambiguously 
represent a covariant spin-3/2 field. The positive energy spinor is defined by the following 
equations: 

(i7-P + M)ii^(P;r) = 0, 7„u^(P;r) = 0, P„it^(P; r) = , (A.IO) 
where r = —3/2, —1/2, 1/2, 3/2. It is normalised: 

u^,(P;r')u,,{P-r) = 2M , (A.ll) 

and satisfies a completeness relation 

3/2 

— MP;r)u4P;r)^A+iP)R,,, (A.12) 

r=-3/2 

where 

1 2 - - 1 ' 

P^f = Sfj.vlD - -7^7^ + -^PiJ.Pi^Io - i-lPfj.^i' - P^7n] I (A. 13) 

with P^ = —1, which is very useful in simplifying the positive energy A's Faddeev equation. 
Appendix B: Nucleon-Photon Vertex 

In order to explicate this vertex, we write the scalar and axial- vector components of the nucleons' 
Faddeev amplitudes in the form [cf. Eq. 1151 1 



^(k-P) 



■ W^'ik-P) ' 




. %{k;P) . 





5(fc;P)u(P) 
A'^{k-P)u{P) 



i = l,...,4. (B.l) 



For explicit calculations, we work in the Breit frame: P^ = P^^ — (5n/2, P^ = P^^ + Qm/2 
and P^^ — (0,0,0, i\/ A/2 + Q2/4), and write the electromagnetic current matrix element as 
[cf. Eq. 1331] 



P'\J""\P) = yl+(P') 

" BF 



-^^Ge + M„-^{Ge - Gm) 



yl+(P) (B.2) 



/ (B.3) 
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In Fig.Qwe have broken the current, J^™ (p, P'; k, P), into a sum of six parts, each of which we 
now make precise. 

The one-body term of Diagram 1 is expressed as 

jr = S(p,)rr (p,; k,)S{k,) (ks) + A'* (fc,)) {27t)^5\p - k - f,Q) , (B.4) 

where r^'^{pq; kg) — Qq rf^{pq; kq), with Qq = diag[2/3, —1/3] being the quark electric charge 
matrix, and rfj,{pq \ kq) is given in Eq. I48II . 
Diagrams 2 and 4 are obtained through 

4" = A'(pi)[r^^''{pa;ka)Y A'{ka)S{kq){2^f5\p-k + riQ). (B.5) 

For Diagram 2: [r^«(pd; fc^)]*^ = diag[Oo+ , Qi+^^m^I. where Qo+ = 1/3 and T/?^ is 
given in Eq. (|^, and Qi+ = diag[4/3, 1/3, — 2/3] with F^^ given in Eq. ||^; while for Di- 
agram 4: [r;?''(pd;fcd)]'^^ = 0, and [/"^"(pd; fed)]^'^ = Tsa, which is given in Eq. JSIJ, and 
[/^'(pd; fcd)]^'^ = -Tas- Naturally, the diquark propagators match the line to which they are 
attached. 

In Eqs. 1B.4II & IB.St the momenta are 

kq^-qP + k, Pq = riP' +p, 

kd^fjP-k, pd = -nP'~p, ^ 

with rj + fj — 1. The results reported herein were obtained with rj = 1/3, which provides a 
single quark with one-third of the baryon's total momentum, but, as our approach is manifestly 
Poincare covariant, the precise value is immaterial. Nevertheless, numerical results converge 
more quickly with this natural choice. 

The quark exchange contribution to the vertex. Diagram 3, is obtained with 

= -lsikq)A\k,)r\p^,ka)S^iq)r^r^iq',q)S^iq')r'^{p'„p,)A'{p,)SiPq) , (B.7) 

wherein the vertex 7^" appeared in Eq. IB. 41 . The full contribution is obtained by summing over 
the superscripts which can each take the values 0"*", l"*". The so-called seagull contributions. 
Diagrams 5 & 6, are given by 

= ^S{kq)A^k,)[x;{pq,q',k,)S^{q')r'^{p'„p,) 

- rip,,k,)S^{q)Xf,{-kq,-q,pi)) A^pi)S{pq), (B.8) 

where, again, the superscripts are summed. In Eqs. 1B.7II & IB.SII the momenta are 
q — fjP — rjP' — p — k , q = f]P' — r\P ^ P ~ k , 

Pi = {p,~q)/2, p'-2 = {-~kq + q')/2, (B.9) 

pi = (P, - g')/2 , p2^{-kq + q)/2. 

Diagrams 1, 2 and 4 are one-loop integrals, which we evaluated by Gaufiian quadrature. 
The remainder. Diagrams 3, 5 and 6, are two-loop integrals, for whose evaluation we employed 
Monte-Carlo methods. 

References 

1. H.y. Gao, Int. J. Mod. Phys. E 12, 1 (2003) [Erratum-ibid. E 12, 567 (2003)]. 

2. V.D. Burkert and T.-S.H. Lee, "Electromagnetic meson production in the nucleon reso- 
nance region," nucI-ex/0407020 

3. M.K. Jones et al. [JLab Hall A Collaboration], Phys. Rev. Lett. 84, 1398 (2000). 

4. O. Gayou, et al, Phys. Rev. C 64, 038202 (2001). 

5. O. Gayou, et al. [JLab Hall A Collaboration], Phys. Rev. Lett. 88, 092301 (2002). 



30 On Nucleon Electromagnetic Form Factors 

6. J. Arrington, Phys. Rev. C 69, 022201 (2004). 

7. I. A. Qattan et at, "Precision Rosenbluth measurement of the proton elastic form factors," 
|nucl-ex/0410010| 

8. G. A. Miller and M. R. Frank, Phys. Rev. C65, 065205 (2002). 

9. S. Boffi, L.Y. Glozman, W. H. Klink, W. Plessas, M. Radici and R.F. Wagenbrunn, Eur. 
Phys. J. A 14, 17 (2002). 

10. M. G. Fuda and H. Alharbi, Phys. Rev. C 68, 064002 (2003). 

11. S.J. Brodsky, J.R. Hiller, D. S. Hwang and V. A. Karmanov, Phys. Rev. D 69, 076001 
(2004). 

12. B. Julia-Diaz, D.O. Riska and F. Coester, Phys. Rev. C 69, 035212 (2004). 

13. R.T. CahiU, C.D. Roberts and J. Praschifka, Austral. J. Phys. 42, 129 (1989). 

14. H. Reinhardt, Phys. Lett. B 244, 316 (1990). 

15. R.T. CahiU, C.D. Roberts and J. Praschifka, Phys. Rev. D 36, 2804 (1987). 

16. M.S. Bhagwat, A. H6U, A. Krassnigg, C.D. Roberts and P. C. Tandy, Phys. Rev. C 70, 
035205 (2004). 

17. C.J. Burden, R.T. CahiU and J. Praschi&a, Austral. J. Phys. 42, 147 (1989). 

18. H. Asami, N. Ishii, W. Bentz and K. Yazaki, Phys. Rev. C 51, 3388 (1995). 

19. M. Oettel, G. Hellstern, R. Alkofer and H. Reinhardt, Phys. Rev. C 58, 2459 (1998). 

20. M.B. Hecht, M. Oettel, C.D. Roberts, S.M. Schmidt, P. C. Tandy and A. W. Thomas, 
Phys. Rev. C 65, 055204 (2002). 

21. A.H. Rezaeian, N. R. Walet and M. C. Birse, "Baryon structure in a quark-confining non- 
local NJL model," |hep-ph/0408233 

22. C.J. Burden, L. Qian, C.D. Roberts, P. C. Tandy and M.J. Thomson, Phys. Rev. C 55, 
2649 (1997). 

23. P. Maris, Few Body Syst. 32, 41 (2002). 

24. A.W. Thomas, Adv. Nucl. Phys. 13 (1984) 1. 

25. R. D. Young, D.B. Leinweber, A.W. Thomas and S.V. Wright, "Systematic correction of 
the chiral properties of quenched QCD," hep-lat/0111041 

26. E.J. Hackett- Jones, D.B. Leinweber and A.W. Thomas, Phys. Lett. B 489, 143 (2000); 
ibid. 494, 89 (2000). 

27. D.B. Leinweber, A.W. Thomas and R.D. Young, Phys. Rev. Lett. 86, 5011 (2001). 

28. P. Maris and C. D. Roberts, Int. J. Mod. Phys. E 12, 297 (2003). 

29. K. D. Lane, Phys. Rev. D 10, 2605 (1974). 

30. H.D. Politzer, Nucl. Phys. B 117, 397 (1976). 

31. C.D. Roberts and A. G. Williams, Prog. Part. Nucl. Phys. 33, 477 (1994). 

32. J. B. Zhang, P. O. Bowman, R. J. Goad, U. M. HeUer, D. B. Leinweber and A. G. Williams, 
"Quark propagator in Landau and Laplacian gauges with overlap fermions," 

[hep- lat /04 1 0045 j and references therein. 

33. C. S. Fischer and R. Alkofer, Phys. Rev. D 67, 094020 (2003). 

34. M.S. Bhagwat, M.A. Pichowsky, C.D. Roberts and P.C. Tandy, Phys. Rev. C 68, 015203 
(2003). 

35. R. Alkofer, W. Detmold, C. S. Fischer and P. Maris, Phys. Rev. D 70, 014014 (2004). 

36. M. S. Bhagwat and P. C. Tandy, Phys. Rev. D 70, 094039 (2004). 

37. M.A. Ivanov, Yu.L. Kalinovsky and C.D. Roberts, Phys. Rev. D 60, 034018 (1999). 

38. A. H611, A. Krassnigg, C. D. Roberts and S. V. Wright, "On the complexion of pseudoscalar 
mesons," nucl-th/0411065 



R. Alkofer, et al. 



31 



39. CD. Roberts, Nucl. Phys. A 605, 475 (1996). 

40. CD. Roberts and S. M. Schmidt, Prog. Part. Nucl. Phys. 45, SI (2000). 

41. R. Alkofer and L. von Smekal, Phys. Rept. 353, 281 (2001). 

42. C J. Burden, CD. Roberts and M.J. Thomson, Phys. Lett. B371, 163 (1996). 

43. P. Maris and CD. Roberts, Phys. Rev. C 56, 3369 (1997). 

44. S. Capstick and W. Roberts, Prog. Part. Nucl. Phys. 45, S241 (2000). 

45. P. Maris, CD. Roberts and P. C. Tandy, Phys. Lett. B 420, 267 (1998). 

46. H.J. Munczek, Phys. Lett. B 175, 215 (1986). 

47. C.J. Burden, CD. Roberts and A. G. Williams, Phys. Lett. B 285, 347 (1992). 

48. M. Stingl, Phys. Rev. D 34, 3863 (1986) [Erratum-ibid. D 36, 651 (1987)]. 

49. G. Krein, CD. Roberts and A. G. Williams, Int. J. Mod. Phys. A 7, 5607 (1992). 

50. S. Ahlig, R. Alkofer, C. Fischer, M. Oettel, H. Reinhardt and H. Weigel, Phys. Rev. D 64, 
014004 (2001). 

51. M. Hess, F. Karsch, E. Laermann and L Wetzorke, Phys. Rev. D 58, 111502 (1998). 

52. P. Maris, "Electromagnetic properties of diquarks," to appear in Few Body Systems, 
|nucl-th /0409008 

53. P. Maris and P. C. Tandy, Phys. Rev. C 60, 055214 (1999); ibid. C 62, 055204 (2000); 
C.-R Ji and P. Maris, Phys. Rev. D 64, 014032 (2001). 

54. P. Maris and P. C Tandy, Phys. Rev. C 61, 045202 (2000); 

55. M. Oettel, L. Von Smekal and R. Alkofer, Comput. Phys. Commun. 144, 63 (2002). 

56. B. C. Pearce and I. R. Afnan, Phys. Rev. C 34, 991 (1986). 

57. N. Ishii, Phys. Lett. B431, 1 (1998). 

58. M. Oettel, M. A. Pichowsky and L. von Smekal, Eur. Phys. J. A 8, 251 (2000). 

59. J. CR. Bloch, A. Krassnigg and CD. Roberts, Few Body Syst. 33, 219 (2003). 

60. J. S. BaU and T.-W. Chiu, Phys. Rev. D 22, 2542 (1980). 

61. R. Alkofer, A. Bender and CD. Roberts, Int. J. Mod. Phys. A 10, 3319 (1995). 

62. F.T. Hawes and M. A. Pichowsky, Phys. Rev. C 59, 1743 (1999). 

63. A. Salam and R. Delbourgo, Phys. Rev. 135, B1398 (1964). 

64. S.J. Brodsky and J.R. Hiller, Phys. Rev. D 46, 2141 (1992). 

65. P. Maris and P.O. Tandy, Phys. Rev. C 65, 045211 (2002). 

66. M. Oettel, R. Alkofer and L. von Smekal, Eur. Phys. J. A 8, 553 (2000). 

67. S.R. Cotanch and P. Maris, Phys. Rev. D 66, 116010 (2002); ibid. D 68, 036006 (2003). 

68. P. MergeU, U. G. Meissner and D. Drechsel, Nucl. Phys. A 596, 367 (1996). 

69. J. Friedrich and T. Walcher, Eur. Phys. J. A 17, 607 (2003). 

70. S. Galster, H. Klein, J. Moritz, K. H. Schmidt, D. Wegener and J. Bleckwenn, Nucl. Phys. 
B32, 221 (1971). 

71. B. Kubis and U. G. Meissner, Nucl. Phys. A 679, 698 (2001). 

72. T. Fuchs, J. Gegelia and S. Scherer, Eur. Phys. J. A 19, 35 (2004). 

73. J. D. Ashley, D. B. Leinweber, A. W. Thomas and R. D. Young, Eur. Phys. J. A 19, 9 (2004). 

74. R. C. Walker et al, Phys. Rev. D 49, 5671 (1994). 

75. T. Sato and T.-S.H. Lee, Phys. Rev. C 63, 055201 (2001). 

76. J. P. Ralston and P. Jain, Phys. Rev. D 69, 053008 (2004). 

77. M. Oettel, "Baryons as relativistic bound states of quark and diquark," PhD Thesis, Uni- 
versity of Tiibingen, nucLth/0012067 

78. P. Maris and C. D. Roberts, Phys. Rev. C 58, 3659 (1998). 

79. A. V. Belitsky, X.-d. Ji and F. Yuan, Phys. Rev. Lett. 91, 092003 (2003). 



